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THEORY OF HEAT TRANSFER AND HYDRAULIC RESISTANCE 

OF OIL RADIATORS* 
By N, B.. Mariamov 

In the present report the coefficients of heat trans- 
fer and hydraulic resistance are theoretically obtained 
for the case of laminar flow of a heated viscous liquid 
in a narrow rectangular channel. The results obtained 
are applied to the computation of oil radiators, which 
to a first approximation may be considered as made up of 
a system of such channels. In conclusion, a comparison 
is given of the theoretical with the experimental results 
obtained from tests on airplane oil radiators. 



NOTATION EMPLOYED 



v mean velocity 0 f flow of liquid- 

V velocity of flow at any point 

t Q temperature of the wall 

t mean temperature of the liquid at any cross section 

t temperature of the liquid at any point 

tj temperature of the liquid at inlet to channel or radiator 

8 = t - t Q mean temperature difference at any cross section 

9 = t - t 0 temperature difference at any point 

^o = *i ~ * 0 temperature difference at inlet to channel or 
.'■ radiator ' - 

a - heat conduction coefficient - 
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\ coefficient of heat conductivity of liquid 

|JL coefficient of Viscosity of liquid 

h height (width) of channel (passage) through which 
liquid f lows 

djj hydraulic diameter of passage 

L length of channel through which liquid flows 

Uu Nusselt number 

Pe^ Feclet number referred to hydraulic diameter 

Re^ Reynolds number referred to hydraulic diameter 

X. coefficient of hydraulic resistance 



INTRODUCTION 



The problem of oil cooling of airplane engines is a 
very important one at the present time. Only a compara- 
tively short time ago this was not the case because of the 
relatively small quantities of heat required to be dissi- 
pated by the radiator which consequently could be of small 
dimensions. Engine cooling by liquids with high boiling 
points considerably increases the heat transfer in the oil 
(1.5 to 2 times) a fact which leads to an increase in the 
required radiator dimensions. At the same time the aero- 
dynamic fineness of airplanes has reached a stage where 
the air-oil radiator, which projects above the surface of 
the airplane, considerably impairs its aerodynamic effi- 
ciency. A study of the operation of oil radiators, with 
the object of obtaining, computation formulas required for 
a rational design of the radiator under given conditions 
(speed of. airplane, temperature of air and oil, power of 
engine, etc.), thus assumes considerable importance. 

The main object of the present paper is that of deter- 
mining the coefficients of heat transfer of oil radiators 
and their hydraulic resistance to the flow of the oil. 
Since the coefficient of heat transfer from the wall to 
the air is sufficiently well known (references 1,2) the 
problem of determining the heat transfer coefficients of 
oil radiators reduces to that of determining the heat 
transfer coefficients from the oil to the wall. 
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The large viscosity of aviation: oils on the one hand, 
and the small width of the _o.il— flow, passages and small, 
velocity of flow on the other .hand , -permit us. to consider 
the flow of the oil within tho radiator as laminar. 
Under this assumption, we have attempted to determine the- 
oretically the heat transfer coefficient from the oil to 
the wa 11. • 

The problem of heat transfer of a straight round tube, 
through which f lows . a heat ed viscous liquid, was first 
.solved by G-raetz and later by a ntimber of other investi- 
gators. (See references.) Notwithstanding the rather 
rough assumptions which were generally made in the solu- 
tion of this problem, the resiilts for small values ..of the 
Reynolds number, as Kraussold has shown (fig. 1), are in 
sufficiently good agreement- with experiment. The applica- 
tion, however , of these results to the computation of the 
heat-transfer coefficients of oil radiators may lead to 
large errors, for it is not known in what manner the data 
obtained for the round tube may be generalized to tho 
narrow rectangular channels, such • as provided by the pas- 
sages for the liquid in airplane radiators with hexagonal 
tubes. For this reason, following the method of Nusselt, 
we have presented the solution of the problem for a nar- 
row rectangular channel so that the solution may be ap- 
plicable to a radiator (fig. 2) which, to a first approx- 
imation, may be considered as a system of such channels. 

To simplify the solution we have made a number of 
assumptions not entirely corresponding to the actual pro- 
cess of heat transfer in the radiator. As will be shown 
below, however , we obtain in- the given case, as also in 
the case of the round tube, a satisfactory agreement of 
the analytical solution with test resvilts. 

The obtained solution .evidently may be applied to 
radiators with any viscoiis fluid as co olant s, provided the 
flow is laminar. 



2. HEAT TRANSFER FOR LAMINAR FLOW OF VISCOUS FLUID 
" IN NARROW RECTANGULAR CHANNEL 



Let us consider the flow of a heated viscous fluid in 
a straight channel having for its cross section a rectan- 
gle whose base is large relative to its height (fig. 3). 
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We shall choose the position of the coordinate axes as 
indicated on the figure. The motion of the liquid may 
he considered as two-dimensional since the side walls, 
in view of their distance from the x-axis, have only a 
slight effect on the velocity field. 

The problem is that of determining the manner in 
which the coefficient of the heat transfer from the liquid 
to the walls of the channel varies with the velocity of 
flow and with the geometric parameters of the channel 
(h,l). For the solution of our problem, we make the fol- 
lowing assumptions: 

1) The temperature of the channel walls is constant. 

2) The channel wall surfaces are absolutely smooth. 

3) The physical constants of the liquid at any cross 

section are constant. 

4) The temperature of the liquid at the inlet section 

of the channel is the same at all points. 

5) The laminar flow of the liquid is fully developed 

over the entire part of the channel under 
consideration. 

Since laminar flow of the liquid is assumed, the phe- 
nomenon of heat transfer from the liquid to the wall will 
be determined chiefly by the conduction of the heat in a 
direction normal to the motion of the fluid. With the co- 
ordinate axes chosen as in figure 3, we write down the 
fundamental differential equation of the steady temperature 
field 



<: 



d 2 * , af* t 9 1£\ - st at at 
ax 2 ay 3 a z 2 / " ax x ay y az 



where a heat conduction coefficient 

t temperature of the fluid at any point 

x,y,z coordinates of any point 

v x rV y' v z components of the fluid velocity at any point 

Denoting the constant temperature of the channel walls 
by t^ and the temperature difference at a given point 
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by 9 - t - t 0 we obtain from (1) the differential equa-v 
tion of steady heat exchange „ 

W s By 3 Sz 2 / 3x * 3y 3b 

Since the laminar flow of the fluid is in the direction 
of the x-azis, we have, obviously 

v = v =0 
y z 

Further, assuming ?1 r const, that is, that the tempera- 

3 x 

ture gradient of the temperature difference of the fluid 
varies only slightly in the direction of the x-axis, we 

2 2 

may neglect the term "by comparison with the term 2..JL 

3 3x 2 dy 2 
Finally, we have '2.-2.=: 0. The differential equation of 

a z 8 

heat exchange under our conditions thus becomes 



3y 3 x 

From elementary hydraulic considerations, it can easily 
be shown that for laminar flow in a channel the cross 
section of which is a rectangle with base large relative 
to its height, the velocity distribution follows the par- 
abolic law as in the case of the round tube and is expressed 
by the formula 



where v is the mean velocity of the fluid 

and h the height of the channel. 

Substituting in (3) the value obtained for v x and 
setting h/2 = b we shall have 

»!i . & »/ i - y*A 99 



3y 2 \ b /3x 
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(The obtained differential equation determines the temper- 
ature difference as a function of the coordinates; the 
other magnitudes - physical constants of the fluid, rel- 
ative length of channel (l/h) , and mean velocity of flow - 
appearing as parameters. 

Following the method of Fourier, we represent the 
general integral of equation (5) in the form of a product 
of two functions: 

9 = cp (x) ii (y) (6) 

one of which is a function of x alone and the other of 
y alone. Substituting in equation (5), we obtain: 



or 




In the above expression, the left side is a function 
of y alone and the right side of x alone. Since the 
equality is true for all values of x and y it follows 
that the right and left sides of the equation must be 
equal to the same constant magnitude. Denoting the latter 

a 

by - £— the differential equation (5) breaks up into a 
b 2 

pair of ordinary differential equations 




The integral of equation (8) is- 

- - 1 4 * 

cp = De 3 v ■ (11) 
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where D and c are constants determined by the 
following boundary conditions: 

6 (0, y) = 9 Q condition at inlet section of channel 

8 (x, b) = 0 condition at surface of channel 

By a change in the variable y/b = | equation (9) is 
transformed into 

^| + c 2 (l- |V = 0 < 12 > 

The solution of the above equation is a certain function 
\}/ (c a ;|)« By applying the method of successive approxima- 
tions to equation (12), we can find this function. Assum- 
ing it as known, we obtain the general integral of equation 
(5) in the form 

_ C 2 2 a x 

9 = De * 3 ▼ \J/(c S , 0 (13) 

The constants of the integration c s and D are deter- 
mined from a consideration of the boundary conditions. 
The constant c s is easily found from the second condition 
which, evidently, can be satisfied only if 

\t/(c S ; 1) = 0 (14) 

In this way we arrived at the equation by which the con- 
stant c 2 is determined. 

The function ^(c 2 ;|), appearing as the solution of 
equation (12), cannot be represented in finite form and 
is expressed in the form of an infinite series. Equation 
(14) thus has an infinite multiplicity of roots. By ap- 
plying the aiethod of successive approximations, we com- 
puted the first three roots 



\|/ (c 2 , 1) x = c j 2 = 2.827 
\|/ (c 2 , l) s = c s 3 = S2.0 
\J/ (c 2 , 1) 3 = c 3 2 = 92.1 
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Corresponding to these, expression (13) for the temperature 
difference may "be represented in the form of an infinite sum 
of particular solutions 

Remembering that 

t>— 2 . « & 
and introducing the Pec.let number 



a 



we obtain 



The constants Dj_ are determined from a consideration of 
the conditions at the inlet section of the channel. For the 
first three constants the following values were obtained 

= 1,208 8 0> 
Do = 0,299 6 0 , 
D s = 0,118 8 9 . 

Substituting the values thus found for c^ 2 and we 
finally obtain 

O.+.H*/''"^. (£+0.299 / 85 '**\ + 

(i) + ■•■]■ <>*> 



- 246 P5T 



+ 0,118 <? t 

On figure 4 are plotted the graphs of the functions \J/ T »^ 8 » 
and ty 3 . To determine the coefficients of heat transfer of 
the channel, it is necessary to find the value of the mean 
temperature difference at the channel section under consid- 
eration. By definition, we have: 

jtivdF 

- ■ *F 

e = 



jvdF 

F 

where v is the velocity of the flow at the point considered 
9" the mean temperature difference at the section considered 



, A 

In our ease +y 



/ 9 (' 



T 2 



dy 
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or i 

/ 6(1 



Since 



therefore 



o 

Suhst itut ing the value of Q from (15), we ohtain: 

1 X 1 

e = e 0 (i,208e 1Mji ~^. J (i-^ i(S)d j=_ 

o 

i * 1 

-85,4 -i" 4" Q /• 

- 0,299 e f/d (5) + 



-246-14- 1 



+o,ii8« PeA 4"j" O-^fcG)*--..). 

n / 



On carrying out all the indicated operations, we ootain the 
following expression for the mean temperature difference 



/ _7 54_L.£- _854 i— -246 — — \ 

f=e 0 (0,916« ' Peh + 0,053 e ' Pe *+ 0,0086 <? Pe " + ...). (16) 

We can now find the expression for the coefficient of heat 
transfer from the liquid to the channel wall. By definition 



where <4 is the heat flow in a direction normal to the 
channel walls. But, as is known, 

where A is the coefficient of heat conductivity of the 
liquid. Hence 



^ 07) 

The derivative 99/9y is obtained from equation (15) as 
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+ 2,242* ' + 1,335, Pe * + ...;. (18) 

Substituting this value in (l?) and the value of 9 
from (16), we obtain 

_ 7 ,54.Li -w^i -2«sT 

X 3,444 e + 2,242 e +1,335. +■• . 

B== "* -7,54 i-f- ' =^4? = ^il- ' 

0,916 e Pe " + 0.053e " " + 0,0086 . *+... 
Dividing numerator and denominator of the ahove expression 

-7.54 - — — 

•by e p * , *e obtain 

X 3,444+ 2,242- + l,335e +... < f2m 

cc = — : — • i 3 I z y ' 

h _ 92,94 i-f- -253,54^ 

0,916+0,053* +0,0086 <? +... 

For comparison of the obtained solution with that of 
Kusselt for the round tube f we introduce the hydraulic 
diameter defined by the formula 

where f 1 is the cross-sectional area of the channel 

a the perimeter of the channel 

Tor a narrow rectangular channel whose base B is large 
in comparison with its height h, the hydraulic diameter 
is equal to 

since the value of h/B may be neglected in comparison 
with unity. 

Introducing the Nusselt number 
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where Fe^ is the peclet number computed for the hydraulic 
diameter d^ = 2h. 

For a round tube, the analytic solution is of the form 

-103 .6 JL Z. -as e JL 5, 

Fe d Pe d 

2.996 + 2.228 e + 1.006 e + ... 



ffu = 



-103 ,8 JL i -336. 6_1_ ■£ 

0. 819 + 0.0976 e Pe. d + 0.01896 e Pe d + . . . 

(See, for example, reference 3, pp. 225-228.) (22) 

Figure 5 gives the curves of the Efusselt number as a 

function of the reduced length -i- 5. for the channel and 

Fe d 

tube. as shown by these curves for both the channel and 
the tube, the intensity of the heat exchange is infinitely 
large at the initial sections. The coefficient of heat 
transfer then decreases in the flow direction very rapidly 
at first, then at a slower rate,, gradually approaching a 
certain limiting value. Denoting the latter by Nu min , 
we obtain from formulas (21) and (22), setting 

1— - = CO 

Pe d 

for the round tube: Nu„. = S.65 . 

min 

for the narrow rectangular channel: 

Nu mia = 7 ' 54 
or the corresponding values of ^min* 

a , = 3. 65 
min d 

a s >'8.77 ^ » 
min h 

He have thus obtained the result that, the minimum 
values of the heat- transfer coefficients for the round 
tube and the narrow rectangular channel are very close to 
each other if the diameter d of the tube is equal to the 
height h of the channel. The length of the "initial part," 
along which the minimum value of the coefficient of heat 
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transfer from the liquid to the wall is practically estab- 
lished, may he found from formulas (21) and (22). If the 
value of the heat- transfer coefficient is permitted to 
differ by 1 percent from the minimum, we obtain for the 
round tube 

| = 0.05. Pe 
for the narrow rectangular channel 

f~ = 0.0125 Pe h 
d h n 

Introducing the parameter 

Pe = t* 
a. 

and remembering that d^ = 2h, we obtain 

~ = 0.05 Pe (23) 

that is, the length of the "initial part" of the narrow 
rectangular channel is equal to the length Of the initial 
part of the round tube if the characteristic linear dimen- 
sions; namely, the height h and the tube diameter d 
are equal to each other. 



3. ANALYSIS OF THE OBTAINED SOLUTION 



The analytical solution was obtained on the basis of 
a number of simplifying assumptions. The latter, as was 
remarked above, do not entirely correspond to the strict 
conditions of the process and the inaccuracies introduced 
in the course of obtaining the solution evidently show up 
in the final result. We shall consider a little more in 
detail the assumptions made and investigate their effect 
on the .final result. An essential simplification of the 
problem is attained by the assumption that d 2 9/dx s is 
negligibly small by comparison with 3 2 9/dy s . As, however, 
may be seen from the formulas obtained, the temperature 
near the initial cross section of the channel changes very 
■rapidly along the x-axis and hence the axial gradients 
here are very large. The assumption, therefore, does not 
correspond to the true character of the heat exchange in 
the initial part* For the elements sufficiently removed 
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from the initial section, the axial gradients are very- 
small and negligible and do not introduce any considerable 
error in the final result. The error is thus very small 
for the region of stabilized heat exchange but may be very 
large in the initial part. It is evident that the smaller 
the initial part as compared with the entire length of the 
channel the smaller will be the error. Since the length 
of the initial part is proportional to Pe the error will 
decrease with decrease in Fe and with increase in the 
relative depth of the channel x/h. 

A second very important assumption is that of the non- 
variance of the physical constants of the viscous liquid 
over the channel cross section. This means that the vis- 
cosity of the liquid is constant over every section and 
that the velocity distribution over any section follows 
the parabolic law. Actually the viscosity of the liquid 
varies considerably with the temperature - the viscosity 
of mineral oils, for example, varying almost in inverse 
proportion to the cube of the temperature. Consequently, 
the actual velocity distribution over a cross section will 
not follow the parabolic law. On figure 6, taken from 
the work of Greber and Erk (reference 3), are shown the 
velocity distributions for laminar flow in a tube for 
three cases: isothermal flow (curve I), case of heat ab- 
sorption by the liquid (curve III), case of heat emission 
by the liquid (curve II). As may be seen from the figure, 
the deviation introduced by the heat exchange in the ve- 
locity distribution from that corresponding to isothermal 
laminar flow is very large. For this reason, the obtained 
theoretical formulas are applicable only if the temperature 
differences are small or, if this is not the case, only to 
liquids the physical properties of which do not change 
with the temperature. 

With regard to the assumption of constancy of the 
temperature of the channel walls, cases may occur in prac- 
tice where this assumption actually applies. Such cases, 
however, are very rare. Generally the temperature of the 
walls changes along the direction of flow of the liquid. 
The assumption of constant wall temperature then introduces 
some errors in the final result. It is readily seen that 
these errors decrease with increase in the velocity of 
flow of the liquid, that is, with increase in Pe and de- 
crease in the L/h, the relative depth of the channel. 
The result is thus obtained that the errors arising from" 
the assumption of constancy of the wall temperature to a 
certain degree are compensated by the errors arising from 
neglecting d s 6/3x 2 . 
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The remaining simplifying assumptions sufficiently 
well correspond to the actual conditions of the heat trans- 
fer and do not greatly affect the final result. 

4. RESISTANCE OF A HEATED VISCOUS LIQUID FOR 
LAMINAR FLOW IN A NARROW RECTANGULAR CHANNEL 

The resistance coefficient for isothermal flow of a 
viscous liquid in a narrow rectangular channel is expressed, 
as is known, "by the following formula 

X fr = 26 (24) 
" Re 

where Re is the Reynolds number computed by the formula 

He = — --- (25) 
p. 

in which the density and the coefficient of viscosity of 
the liquid are, respectively, denoted by p and |A. 

We shall generalize formula (24) to the case of the 
flow of a viscous liquid whose temperature and, hence, also 
whose physical constants vary over the cross section of 
and along the channel. For this purpose we choose the co- 
ordinate axes as shown in figure 7 and consider an elemen- 
tary parallelepiped of liquid of volume 2ydx 1. Denoting 
by t the tangential stress on the surface of the element, 
we have, from the condition of equilibrium, 

2 ydp = 2Tdx 
dv 

where t = u — 

Remembering that the derivative dv/dy is negative, we 
obtain 

dp dv 
dx dy 

In the case of noni sothermal flow of the liquid the 
coefficient of viscosity \i is a function of x and y. 
Denoting this function by y,(x,y), we substitute it in the 
obtained equation: 

y j£= -n(x.jr) j£ 

dx dy 
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Separating the variables and Integrating with respect 
to y between the limits y and b, we shall have: 

b b 

f dv = - *Z f 
. J dx.J n(x,y) 

y y 

and since v = 0 for y = b 



o 



■f^-T (2 6) 



y 



In order to compute the obtained expression, it is 
necessary to know the temperature field of the liquid at 
the section considered and the law of variation of the 
viscosity with the temperature. The temperature field at 
any cross section of the channel may be determined by for- 
mula (15). The application, however, of this formula to 
the computation of expression (36) is very inconvenient. 
For our purpose it is sufficient to make use only of the 
first term of this formula, giving the temperature field 
in the region of stabilized heat interchange. .In fact, 
as may be seen from figure 8, starting with the value 

1 I 

P~e d = 0,005 t the first term, with sufficient accuracy, 

characterizes the temperature field, the length of the 
initial portion, as follows from formula (23), being equal 
to 

|- = 0.005 Pe. 



In present-day airplane oil radiators, the initial 
portion constitutes about 15 to 30 percent of the entire 
relative length of the channel passage through which the 
oil flows. Thus the greater part of the channel is in 
the region of stabilized heat interchange and hence in 
making use only of the first term of the formula (15) for 
the computation of expression (36), we do not make a large 
error. 

The law of variation of the viscosity with the tempera- 
ture is generally given by an empirical formula. In the 
range of 50° to 150° the change in viscosity of aviation 
oils with temperature is expressed with sufficient accuracy 
by the following empirical formula (reference 4) 
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where i is a number characteristic of the oil 

t the temperature of the oil in degrees C. 

We determine the temperature of the oil at any point 
of the cross section under consideration by (15), taking 
only the leading term of the series 

t = t Q + 1,208 6 o e~ * 54 *e h yj 



The function y x = ^ 2 (y/b), as may be seen from figure 9, 
may with a sufficient degree of accuracy be set equal to 

$) (28) 

_, 1 X 

Then -7.54 / N 

t = t 0 + 1.208 6 0 e pe h (^l - Z±J 



Substituting the found value of t in formula (27) 
and setting 

a -7. 54 .L 5 

e 0 Pe h 

0 = 1.208 ~- e (29) 
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we obtain 

1000* 1 

*=-Tr 



Eoting that 1000 i/.t* is the viscosity corresponding to 
the temperature of the wall, we may write: 

t*o 



where M-o iB the viscosity corresponding to the temperature 

of the wall. Substituting the obtained value of |j, in 
(26) , we obtain: 



V 

We introduce a new variable 



Then 



!*0 

o r 



-(-*) 
— ■=-£-f/<"+ a ** 



o 

Carrying out the computation of this integral, we obtain- 

— fc«[0-*-)+f«( , -f) , + e, ('-4)'+ 



The mean velocity of flow of the liquid is determined 
from the relation 



b b 



2 j'vdy 

0 

or substituting the value of v from (30) 
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The computation of this integral gives 

5 = "|r (1 + 1,2 c + 0,687 ^ + 0,0682 C3) - (31 > 

ffrom this we have 



3p 0 v 1 

~dx = 6 s ' 1 + 1,2 C+ 0,687 C* + 0,0682 C 3 

Sett ing 

» =1,208 4*-- 

then 



(32) 
(33> 



_ 75 4_L_ Jl 

C = «e ' Pe * (34) 

and 

4p _ 3p 0 y 1 



P _754 — — / -7S4-L JL\ S / _7 54 _ 1 _jfA :i 

1 + l,2ne * + 0,687 \ n<? * * * J +0,0682 [tie P " ; 

integrating this expression between the limits 0 and I, 
we obtain 



/*- 



— & J 



dx 



J _754J_.fl / _7 54_LjL\ S / _ 754 _LA\* 

0 l + l,2n<? ' P ' ft + 0,687 Pe ft J +0,0682 (^ne ' Pe * j 



o r 

■=lPi—Po) = 

_ 3y. 0 v f dx_ 

b 2 J ~ 



J -754 — — / - 754-L-iV / _7 S 4±JL\ 3 " 

0 \ + \,2ne ' "» +0.687 ' * J + 0,0682 (^e ' Pe " ) 



where 



Ap is the resistance of the channel for the flow of 
heated oil 

p 0 pressure at the initial section of the channel 

p^ pressure at the final section of the channel 

To compute the obtained integral we make the substitution 
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754 — ^- 
/,S4 Pe h 



and set 



Then* 



-= 7 ' 54 ^ =30 : 2 ^:i- (35 > 

ne~ m 

. SppPehv r , dz 

Ap== 6*7,54 J z (1 + 1 ,2 2+ 0,687 z* + 0,0682 z«) - 

5 

♦The parameter Fe changes only slightly along the chan- 
nel so that with sufficient accuracy Pe may "be considered 
independent of the length. 

Carrying out the integration, we obtain after certain 
transformations _ 

3h. 0 Pehv [ . nnnA ^,_ ne m +8,15 

0,682ra 2 g- 2/ " + 0,132ftgr m + 0,12 3 
+ °' 51n 0,682 « 2 + 0,132« + 0,123 " r 

2,32 «(g~ m — 1) 1 

• + 0,316 afC tg 1 + (0,23 + 2,32 ne~ m ) (0,23 + 2,32 n) J ' 



By definition, we have 

p 2 2A 



^r = -^-T (36) 



Hence 



fr 7,54 P t»6 2 /L ■ n + 8,15 



,15 

0,682 tt 2 e~ 201 + 0, 1 32 ne~ m + 0, 1 23 , 
' 0,b82rt 3 + 0,132« + 0,123 " ' + 

+ 0,316arctg 2 ' 32 ^"l" *> 1 

■ 1 + (0,23 -f- 2,32 ne~ m ) (0,23 + 2,32 n) J 

Heme'mbering that 

h 2 Peh Pe h a b 1 
* 2=?? "T ; 77547 = 30^27 = 

and denoting 

^ 0 = 2 ^,. : ' (37) 

we obtain finally 

\ 96 f i nnrnvi «g~ m + 8,15 , ■ 

\ = R^t [ « + 0,0047 In ,. + 

, _ 0,682nV~ 2m + 0, 1 32 ne~ " + 0, 1 23 , 
+ U,t>ln- 0 682w 2 + 0,132 ra + 0,123 . i_ 

+ 0,3t6arctg 2,32,^-1) j 

^ 1 + (0,23 + 2,32«g-' n ) (0,23 + 2,32 n) J 
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In the particular case thai; the flow is isothermal, 

e 

96 



we have 9 0 = 0 and therefore n = 0 , whence 



that is, ire obtain formula (24). 

The numerical value of the term 



0 .00 47 In 



ne" m + 8.15 
n + 8.15 



is very small by comparison with the other terms entering 
expression (38) and may be neglected. 

Setting 

. _ _ . 0.68 2 n 3 e"* sm + 0.132 ne"* m + 0,123 
M = m + 0.5 in " — 

0.682 n s + 0.132 n + 0.123 

i* 2.32(e" m - 1) , x 

+ 0.316 (tan -1 ) (39) 

1 + (0.23 + 2.32ne~ m ) (0.23 + 2.32n) 



we obtain K r = f 6 - & (40) 

R 0 m 



where B. Q is the Reynolds number referred to the tempera- 
ture of the wall. 

The formula obtained permits the computation of the 
resistance coefficient for laminar flow of a heated oil in 
radiators with oil passages in the form of narrow rectangu- 
lar sections whose base is large by comparison with the 
height. The formula may, however, be employed also for 
other shapes of pas sages .provided that the height of the 
channel is small compared to the base, since for this con- 
dition the shape of the channel plays an insignificant part. 
For the computation of the hydraulic resistance of a radia- 
tor, it is necessary to know the mean temperature t 0 of 
the radiator walls, the inlet temperature of the oil t^ 0 

and the reduced length of the channel r 5 --.. This 

* e h d h 
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formula" evidently holds only under the condition that the 
temperature of the" oil throughout the channel is not be- 
low 60° to 50° 0 since at lower temperatures the law as- 
sumed by us for the variation of the viscosity "becomes un- 
true and the formula gives too great a value for the co- 
efficient Xf r * It is not difficult to obtain a formula 
for 'application also to lower oi 1. t emp erature s . For this 
purpose it is necessary to assume another law of variation 
of viscosity with temperature: namely,. one that corresponds 
to the new temperature range. In airpiane oil radiators 
the temperature of the oil as a rule is always above 50° C. 
In view of this fact. we considered it sufficient to limit 
ourselves to the derivation of a formula for the computa- 
tion of the resistance for an oil temperature above 50° 0. 
For convenience in .applying formula (40) to practical com- 
putations the chart of figure 10 was constructed from which 
the value of \f r for given values of n and m can 
readily be found. 



5. COMPARISON OF THEORETICAL WITH TEST RESULTS 



In the testing of airplane radiators the coefficient 
of heat transfer from the liquid to the air is generally 
determined from the relation 

k = — _3 ( 41 ) 

S(H - t a ) 

where k is the overfall coefficient of heat transfer in 
cal/hr • 

S the cooling surface of the radiator in m s 

t^ the arithmetic mean of the inlet and outlet tem- 
peratures of the cooling liquid in the radiator 

t a the temperature of the air entering the radiator 

The coefficient of heat transfer from the liquid to 
the wall is separated from the over-all coefficient by 
various methods. A typical method is the following. The 
heat transfer is determined for the same radiator for two 
different liquids (for example water and ethylene glycol 
or water and oil) for the same weight of air through the 
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radiator. Since the coefficient of heat transfer from the 
wall to the air should he the same for both cases the co- 
efficient of heat transfer from the liquid to the. wall can 
be found from a comparison of the tests. 

As may be seen from the method itself of determining 
the heat-transfer coefficient of the radiator from the 
liquid to the wall, a mean coefficient is obtained for the 
heat given off over the entire cooling surface in contact 
with the liquid referred to the temperature difference be- 
tween arithmetic mean of the inlet and outlet temperatures 
and the mean temperature of the radiator cooling surface. 
Hence for the comparison of the obtained theoretical re- 
sults with experimental results we must average the heat 
transfer coefficients for the entire cooling surface of 
the channel and refer to the temperature difference cor- 
responding to the test data. For this purpose we compute 
the quantity of heat which the heated liquid gives off in 
laminar flow in a rectangular channel whose length is equal 
to L and perimeter unity. If ^ is the quantity of 
heat in cal/hr 

ft = J <X 9 dx (42) 
0 

where a is the coefficient of heat transfer from the 
liquid to the wall at the channel cross 
section under consideration. 

0 the mean temperature of the liquid at the section 

Substituting the values of "9 and a from equations (16) 
and (20) and retaining only three terms of the series. we 
obtain 
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which after integrating gives 

Q = w 0 I 0,49 — 0,458 e ' Pe h — 0,0263 e * * — 



- 0,00543 e^^)- < 43 > 



The mean coefficient of heat transfer from the liquid to 
the wall will be equal to 

r_ Q 



(44) 



where 9^ is the mean temperature difference at the end 
of the channel. Prom (16) we. obtain for the expression 



G £ = K (o,916e 7,54 p» * + 0,053e 85,4 ^ * -f 0,0086<? 246 p ' * + . . .;. 



(45) 



Substituting the values obtained for C4 and 9" L in (44) 

and retaining for 9~ L only the first two and for Q only 

the first three terms of the series, as is sufficient for 
our purposes, we obtain after some transformations 

30 2 ,i 341 \ ^ 

7—p„ L k 0,49^ 0,458g ' p ' h " h — 0,0263* Pe *> rf «< 

« r j T j j 

^ -30,2 -J- -i~ -341 — — 

l + 0,916e Pe h "k" + 0,053e p ' h * h 
or introducing the Nusselt number 

/Vw = , 

we obtain 

-30,2 -'- — - 341 1 — 

JV„ = p. A 0,4 9-0,458g P ' h * h — 0,0263g 

^ -30.2 » * -wi-I— L. ' ( 46 ) 

l + 0,916e Pfh rfh 4-0,053e p ' h " h 



Figure 11 gives tee curves of Nu against Pe^ d^/L 

drawn by the above formula. On the same figure is given 
the curve constructed from the t^est results of Brown and 
Barlow (reference 1) for the coefficient of heat transfer 
from ethylene glycol to the walls of a honeycomb radiator. 
The latter was assembled from hexagonal tubes in such a 
manner that at two out of every six sid.es of a tube wide 
passages h were formed and at the remaining four sides 
narrow ones h» where 2h» = h. - Thus the velocity of 
flow of the liquid in any cro ss- sect ion of the radiator 
was practically constant. The hydraulic diameter for 
the liquid passages was computed by the formula 
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6 1+ 2h + 4h« 



(47) 



where I was the length of the air tube (depth) of the 
radiator. 

As may be seen from figure 11, the theoretical curve 
lies sufficiently close to the experimental curve. 

The same good agreement is obtained from a comparison 
of the theoretical results with the test results for air- 
oil radiators. The latter (fig. 13), are honeycomb radi- 
ators with round tubes. To form the passages for the oil, 
the ends of the tubes are formed into a hexagon and soldered 
together. In contrast to the usual honeycomb radiators the 
coolant (oil) flows not in a direction perpendicular to 
that of the air flow in the pipes but parallel to it. 

With the object of attaining a more uniform operation 
of oil radiators, partitions are constructed by means of 
which the velocity of the oil flow is varied at different 
sections of the radiator. To compute the reduced length 
of the radiator and the Nusselt number, it is necessary to 
find the averaging law of the velocity and the expression 
for the hydraulic diameter. Denoting by v and 3P with 
the corresponding subscripts the velocity of the liquid 
and the cro s s- sect ional area of the passages in one sec- 
tion, the mean velocity of the liquid in the entire radia- 
tor will be 



v^ + v 2 



v V 
n n 



where 3P is the over-all area of the cross sections of all 
passages for the liquid. From elementary con- 
siderations we obtain 



where n 
W 



v = 



(n + 1) W 



3.6X10 6 F (1 - f) 
is the number of partitions 
flow of oil through the radiator in m 3 /hr 



(48) 



frontal area of the radiator in 



m 



2 
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f the coefficient of "live" cross section of the 

radiator equal to the ratio of the air-passage 
area to the entire frontal area of the radiator. 

Likewise from elementary considerations we may obtain the 
following expression for the hydraulic diameter of the oil 
passages 

i.xq» (ytjni.., (49) 

h d 0 a ° 

where d Q is the outside diameter of the tube in m 

h^ the minimum oil passage width 

Figure 14 gives the curve of Nu against Pejj d h /l> 
constructed by formula (46). There are also plotted the 
test points for air-oil radiators obtained from the radia- 
tor tests of S. P. Cherbakov at CAHI . The velocity of 
flow of the oil and the hydraulic diameter were computed 
by the above formulas. 

Formula (40) obtained for the resistance of a honey- 
comb radiator with heated-oil flow gives satisfactory 
agreement with test results as may be seen from figure 15, 
which gives the variation of the resistance coefficient \^ 
with the Reynolds number Re 0 . In computing ~\f T by for- 
mula (40) the temperatures taken for the radiator walls 
were those computed from the test data of S. F. Cherbakov 
(reference 9) . 

CONCLUSIONS 



1. The theoretical values obtained for the coefficient 
of heat transfer from the liquid to the wall for laminar 
flow of the liquid in a rectangular channel for Pe^ dj 1 /L<150 
g;jLve satisfactory agreement with the test results. 

2. The theoretical values obtained for the resistance 
coefficient likewise give satisfactory agreement with test 
results for honeycomb radiators in which the flow of the oil 
is directed perpendicular to the direction of the flow of 
the air in the tubes. 

In conclusion the author wishes to express his sincere 
thanks to 0. X. Blasov for suggestions of value in prepar- 
ing this paper. 
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Figs. 1,3,03 




Figure 1.— Theoretical curve of Ru against |pe for round 
tubes for parallel velocity distribution* 




Figure 3.— Arrangement of tubes and channels for passages 
for the liquid in a radiator system. 




Figure 18* 
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Figure 5.««» Curves of Nu against l/Pe-x/d. 
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Figure 6.- Velocity distribution in pipe for laminar flow. I- isothermal 

flow, II- emission of heat by the liquid, III- absorption of 
heat "by the liquid. Mean velocity v in all three cases the same. 
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Figure 7 
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Figure 8.- Curves of 0/0q against y/t>. Continuous curve computed "by- 
formula (15) , dotted curve "by first term of formula. 
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Figure 10.- Chart for computing the coefficients of hydraulic resistance 
of oil radiators. 
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Figure 11.- Comparison of theoretical with test results (Brown-Barlow) 
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Figure 15.- Comparison of theoretical curve of X f r against Re 0 with 
test results. 
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